Molecular insights into the inhibition mechanism of harringtonine against essential proteins associated with SARS-CoV-2 entry.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently posed a serious threat to global public health. Harringtonine (HT), as a small-molecule antagonist, has antiviral activity against a variety of viruses. There is evidence that HT can inhibit the SARS-CoV-2 entry into host cells by blocking the Spike protein and transmembrane protease serine 2 (TMPRSS2). However, the molecular mechanism underlying the inhibition effect of HT is largely elusive. Here, docking and all-atom molecular dynamics simulations were used to investigate the mechanism of HT against the receptor binding domain (RBD) of Spike, TMPRSS2, as well as the complex of RBD and angiotensin-converting enzyme 2 complex (RBD-ACE2). The results reveal that HT binds to all proteins primarily through hydrogen bond and hydrophobic interactions. Binding with HT influences the structural stability and dynamic motility processes of each protein. The interactions of HT with residues N33, H34 and K353 of ACE2, and residue K417 and Y453 of RBD contribute to disrupting the binding affinity between RBD and ACE2, which may hinder the virus entry into host cells. Our research provides molecular insights into the inhibition mechanism of HT against SARS-CoV-2 associated proteins, which will help for the novel antiviral drugs development.

Identification of Potential Lead Compounds Targeting Novel Druggable Cavity of SARS-CoV-2 Spike Trimer by Molecular Dynamics Simulations.

The global pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become an urgent public health problem. Spike (S) protein mediates the fusion between the virus and the host cell membranes, consequently emerging as an important target of drug design. The lack of comparisons of in situ full-length S homotrimer structures in different states hinders understanding the structures and revealing the function, thereby limiting the discovery and development of therapeutic agents. Here, the steady-state structures of the in situ full-length S trimer in closed and open states (Sclosed and Sopen) were modeled with the constraints of density maps, associated with the analysis of the dynamic structural differences. Subsequently, we identified various regions with structure and property differences as potential binding pockets for ligands that promote the formation of inactive trimeric protein complexes. By using virtual screening strategy and a newly defined druggable cavity, five ligands were screened with potential bioactivities. Then molecular dynamic (MD) simulations were performed on apo protein structures and ligand bound complexes to reveal the conformational changes upon ligand binding. Our simulation results revealed that sulforaphane (SFN), which has the best binding affinity, could inhibit the conformational changes of S homotrimer that would occur during the viral membrane fusion. Our results could aid in the understanding of the regulation mechanism of S trimer aggregation and the structure-activity relationship, facilitating the development of potential antiviral agents.

High-Throughput Screening for the Potential Inhibitors of SARS-CoV-2 with Essential Dynamic Behavior.

Global health security has been challenged by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic. Due to the lengthy process of generating vaccinations, it is vital to reposition currently available drugs in order to relieve anti-epidemic tensions and accelerate the development of therapies for Coronavirus Disease 2019 (COVID-19), the public threat caused by SARS-CoV-2. High throughput screening techniques have established their roles in the evaluation of already available medications and the search for novel potential agents with desirable chemical space and more cost-effectiveness. Here, we present the architectural aspects of highthroughput screening for SARS-CoV-2 inhibitors, especially three generations of virtual screening methodologies with structural dynamics: ligand-based screening, receptor-based screening, and machine learning (ML)-based scoring functions (SFs). By outlining the benefits and drawbacks, we hope that researchers will be motivated to adopt these methods in the development of novel anti- SARS-CoV-2 agents.

Recognition between CD147 and cyclophilin A deciphered by accelerated molecular dynamics simulations.

CD147 functions as the receptor of extracellular cyclophilin A (CypA) in various diseases, and CD147-CypA binding ulteriorly underlies the pathological process of various viral infections including HIV-1, SARS, and SARS-CoV-2. Although CyPA has been identified as a key intermediate pro-inflammatory factor, the mechanism by which CD147 cooperates with CypA in the development of the cytokine storm remains largely unknown, and the binding profile of CD147 with CypA remains to be elucidated as well. Here, we prepared three binding models of the CD147-CypA complex, including the active site of CypA severally binding to the groove bound by the Ig1 and Ig2 domains (model-0), P180-G181 (model-1), and P211 (model-2) of CD147, as well as introducing mutations P180A-G181A and P211A individually in each model. All systems were studied using accelerated molecular dynamics simulations and the molecular mechanics generalized Born surface area (MM/GBSA) method. For model-0, CypA bound to the ectodomain of CD147 with the highest binding affinity. Moreover, mutations P180A-G181A of CD147 in model-0 decreased the binding affinity and weakened the dynamic correlation between CD147 and CypA, which resulted in CypA shifting from the initial binding location. Other residue mutations of CD147 did not significantly affect the CD147-CypA binding, as reflected by the energy and structural analyses. Compared with surface plasmon resonance results and nuclear magnetic resonance shift signals, CypA should tend to reciprocally bind to the groove of CD147, and the binding process might be modulated by P180-G181 rather than P211. Besides, residue R201 of CD147 is critical for CD147-CypA binding and needs further experimental verification. These findings further our understanding of the recruitment between CD147 and CypA and its potential role in the development of inflammation and viral infection.

Molecular insights into the binding variance of the SARS-CoV-2 spike with human, cat and dog ACE2 proteins.

SARS-CoV-2 has recently caused an epidemic in humans and poses a huge threat to global public health. As a primary receptor of SARS-CoV-2, angiotensin-converting enzyme 2 (ACE2) exists in different hosts that are in close contact with humans, especially cats and dogs. However, the underlying mechanism of how the spike receptor binding domain (RBD) of SARS-CoV-2 cooperates with human ACE2 (hACE2), cat ACE2 (cACE2) and dog ACE2 (dACE2) and the variation in binding remains largely unsolved. Therefore, we explored the binding behavior of the spike RBD with cACE2, dACE2 and hACE2 via all-atom molecular dynamics simulations. In accordance with the binding free energies and residue interactions, the spike RBD has respective binding specificities with cACE2, dACE2 and hACE2, and the binding affinities decrease in the order of hACE2, cACE2, dACE2, mainly due to changes in the amino acids Q24L, H34Y, and M82T in cACE2 or dACE2. Furthermore, alanine scanning analysis results validated some key residues of the spike RBD interact with ACE2 and provided clues to the variation of amino acid that could influence the transmissibility or immune responses of SARS-CoV-2. Decreasing dynamic correlations strengths of ACE2 with the RBD were found in all hACE2-RBD, cACE2-RBD and dACE2-RBD systems. The ACE2 protein shows variable motion modes across the zinc metallopeptidase domain, which induces different interactions between ACE2 and the RBD. Our studies reveal that the motion pattern of the zinc metallopeptidase domain is critical to the binding behavior of RBD with ACE2. These findings could aid our understanding of selective recognition involving various ACE2 with the SARS-CoV-2 spike and shed further light on the binding mechanisms.